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ABSTRACT. In this paper we study a mathematical model of one-dimensional swimmers per- 
forming a planar motion while fully immersed in a viscous fluid. The swimmers are assumed to 
be of small size, and all inertial effects are neglected. Hydrodynamic interactions are treated 
in a simplified way, using the local drag approximation of resistive force theory. We prove ex- 
istence and uniqueness of the solution of the equations of motion driven by shape changes of 
the swimmer. Moreover, we prove a controllability result showing that given any pair of initial 
and final states, there exists a history of shape changes such that the resulting motion takes 
the swimmer from the initial to the final state. We give a constructive proof, based on the 
composition of elementary maneuvers (straightening and its inverse, rotation, translation), 
each of which represents the solution of an interesting motion planning problem. Finally, we 
prove the existence of solutions for the optimal control problem of finding, among the histo- 
ries of shape changes taking the swimmer from an initial to a final state, the one of minimal 
energetic cost. 
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1. Introduction 

In this paper we study the self-propelled planar motions of a one-dimensional swimmer 
in an infinite viscous three-dimensional fluid. We are interested in the swimming strategies 
of small organisms that achieve self-propulsion by propagating bending waves along their 
slender bodies (such as, for instance, sperm cells and Caenorhabditis elegans). At these 
length scales, viscosity dominates over inertia: accordingly, we ignore all inertial effects in 
our analysis. 

The study of the self-propulsion strategies of microscopic living organisms is attracting 
increasing attention, starting from seminal works by G. I. Taylor II25II . M. J. Lighthill II20II . 
E. M. Purcell [23], and S. Childress |8]. We refer the reader to the recent review lfl"9H for a 
comprehensive list of references. Among the recent mathematical contributions we quote 
Ifl8l [l3l 1241 [5j [2l [4). Many of these papers approach swimming problems within the 
framework of control theory, and this is exploited in (HSU for the numerical computation of 
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energetically optimal strokes. While the connection between swimming and control theory 
is very natural, only recently has this point of view started to emerge and become widely 
appreciated, see IfTUl and the other chapters in the same volume. 

When inertial forces are neglected, and external forces such as gravity are not present 
(neutrally buoyant swimmers), the equation of motion for a swimmer become the state- 
ments that the total viscous force and torque exerted by the surrounding fluid vanish. In 
order to take advantage of the simplifications deriving from the special one-dimensional 
geometry of our swimmers, we adopt here the local drag approximation of Resistive Force 
Theory, first introduced in 03), then also used in E2H and further discussed in IfTTH . It is 
a classical and popular theory widely spread among biological fluid dynamicists, which 
has recently been proved to be accurate and robust in the study of the motion of one- 
dimensional bodies in the length scales and regimes we are interested in, as it is shown, 
e.g., in II 121 . According to resistive force theory, the external fluid exerts on the swimmer 
a viscous force per unit length which, at each point of the swimmer, is proportional to the 
local tangential and normal velocities at that point, through positive resistance coefficients 
denoted by C T and C v , respectively. 

For every t in the time interval [0,T], let s n- x( s >*) De the parametrization of the 
swimmer position with respect to an absolute external reference frame (lab frame), where 
s e [0,L] is the arc length parameter. It is possible to factorize this function as x( s >*) = 
r(t) o£(s, t) , where r(t) is a time dependent rigid motion and s m> £(s, t) describes the shape 
of the swimmer at time t with respect to a reference system moving with the swimmer 
(body frame). 

We suppose that the shape function £ is given. The first problem we address in this paper 
is to determine the rigid motion t ^ r(t) that results from a prescribed time history of 
shape changes t H- £(s, t) . This is obtained by imposing that x = r °£ satisfies the equations 
of motion (the resultant of viscous forces and torques generated by the interaction between 
the swimmer and the fluid vanish for every t ) and solving the resulting force and torque 
balance for r in terms of the given £ . 

Our main result on this first problem is that, if £ satisfies suitable regularity condi- 
tions which are listed in the hypotheses of Theorem 13.31 then the rigid motion r(t) can 
be determined as the unique solution of a system of ordinary differential equations in the 
independent variable t. Therefore, for every initial condition r , there exists a unique r(t) 
such that the resulting function x(s, t) = r(t) o £(s, t) satisfies the force and torque balance. 
In other words, Theorem 13.31 states that looking for a motion that satisfies the force and 
torque balance is equivalent to assigning the shape function and solving the equations of 
motion. 

The second problem we address in this paper is that of controllability. Given a time 
interval [0, T] and arbitrary initial and final states of the swimmer described by the arc 
length parametrizations s n- Xm(s) and s <-> Xfin(s), can we find a self-propelled motion 
x(s, t) in the lab frame such that x(s, 0) = Xin(s) and x(s, T) = Xfin(s) ? By a self-propelled 
motion we mean one such that the equations of motion are satisfied which, in the case of 
self-propulsion, reduce to the vanishing of the total viscous force and torque. The answer 
is affirmative and is contained in Theorem l4.il Our proof is constructive. Indeed, we ex- 
hibit an explicit procedure to transfer Xin onto Xfin based on the composition of elementary 
maneuvers: straightening of a curved configuration and the corresponding inverse maneu- 
ver (i.e., how to map a straight segment onto a given curved configuration), rotation of a 
straight segment around its barycenter, translation of a segment along its axis. Solving the 
motion planning problem for these elementary maneuvers is interesting in its own right, 
independently of the general controllability result, and this is done in Section |4j 

More in detail, given two configurations Xin in Xfin , we show how to straighten them in a 
segment-like configuration, say E in and Sg n , respectively, thanks to Theorem l3.3l Then we 
show how to transfer £ in into £fi n , by explicitly constructing a way to make a rectilinear 
swimmer to translate (without rotating) along its axis, see Section I4TT1 and a way to make 
it rotate (without translating) about its barycenter, see Section 14.21 These constructions 
use suitable bending wave forms that propagate along the body of the swimmer. 

It is interesting to notice, and this will be clear in Section |4l that a very convenient way 
to describe such transformations is by using the angle that the tangent of the swimmer 
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makes with the positive horizontal axis. This angle is given as a function of the time t and 
of the arc length parameter s. This agrees with the traditional approach of prescribing 
the curvature function, since the latter can be recovered by differentiating the angle with 
respect to s (see Remark I3.2D . This classical approach is motivated by the fact that the 
swimmers we are interested in accomplish the shape changes required for force and torque 
balance by relative sliding of filaments along their "spine", hence inducing local curvature 
changes. 

The last problem we address is the existence of an energetically optimal swimming 
strategy. In Theorem 15.11 we prove that, under suitable conditions, there exists a self- 
propelled motion X (s,t) minimizing the power expended. The key hypothesis is a sort 
of non-interpenetration condition for the enlarged body obtained by thickening the curve 
describing the swimmer to a tube of constant thickness. This condition rules out self- 
intersections of the swimmer and yields an a-priori bound on its curvature. 



2. Mathematical statement of the problem 

In this section we describe the mathematical setting for studying the swimming problem 
by adapting to our specific case of a one-dimensional swimmer with a local fluid-swimmer 
interaction the framework introduced and described in H91 I21II . 

Throughout the paper we fix L > to be the length of the swimmer and T > so 
that [0,T] is the time interval in which the motion occurs. We study planar motions in 
three dimensions, and therefore the position of each material point of the swimmer will be 
described by a function x ■ [0, L] x [0, T] — > R 2 , where s e [0, L] is the arc length parameter; 
this request means that for every t the map s n> x (s, t) is Lipschitz continuous from [0, L] to 
K 2 and \ X '(s, t)\ = 1 , where x' '■= dx/ds. As for the derivative with respect to t, x '■= dx/dt 
is intended in the distributional sense as the object that makes the following equality hold 
true 

J J x(MMM) dsdt = - ^ J x (a,t) Mflfl dadf) 

for every tp e C c °°((0,£)x(0,T)). 

We now introduce the local expressions for the line densities f(s, t) and m(s, t) of viscous 
force and torque, as dictated by resistive force theory. Since f(s, t) lies in the plane of the 
motion, m(s, t) is orthogonal to it and is identified with a scalar. They are given by 

f(s, t)= - [C T X r(s, t) X '(s, t) + C v %,(a, t) J X '(s, t)] = -K x (s, t)x(s, t), 
(2.1) m(s, t) = (f(s, t), J X {s, t)) = - (C T X r(s, t) X '(s, t) + C v xAs, i) Jx'(s, t), J X (s, t)) 
= -{K x {s,t)x(s,t),J X (s,t)}. 

Here, C T and C v are positive constants, Xt and x„ are the tangential and normal com- 
ponents of the velocity x, i-e., Xr(s,t) = (x(s,t), x' and Xu(s,t) = (x(s,t), J x '{s,t)) , 



while J 



-1 

1 



is the counter-clockwise rotation matrix of angle tt/2 and 



(2.2) K x (s, t) := C tX '(s, t) ® x '(s, t) + C v {J X '{s, t)) ® (J X '(s, t)) , 

where for any two vectors a, b G M 2 the matrix a ® b is defined by (a ® b)ij = afij . The force 
and torque balance can be written as 



(2.3a) = F(t):= f(s,t)ds = - K x (s,t) x (s,t) da, 

Jo Jo 

(2.3b) = M(t):= m(s,t)ds = - {K x (s,t) X (s,t),J X (s,t))ds, 

Jo Jo 

for a.e. t e [0,T]. 

Remark 2.1. An important remark on the structure of the viscous force and torque is in 
order, leading to a rate independence property. Let tp be a C 1 strictly increasing function 
with ip'(t) > for every t G [0,T]. Then a rescaling in time by tp has no consequences 
on the equations of motion. Indeed, we prove that if x (s, t) satisfies the force and torque 
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balance 02. 3D . then also (x ° y>)(s,t) := x(s, ) does. Let us rewrite the force 02.3aD as 

F x (t) = - J K x {s, t)x(s, t) ds. Then we have 

F xov (t)=- / K xolp (s,(p(t))x(s,<p(t))<p(t)ds 
Jo 

= - <p{t) [ K x (s, <p(t)) X (s, <p{t)) ds = <p(t)F x (cp(t)) = 0, 
Jo 

since F x (tp(t)) = for a.e. t e [0, T] . The same can be obtained for the torque M XO(p (t) . 

We conclude this section by introducing a function space X containing our state func- 
tions, as well as the shape functions: 

(2.4) X:={ X : [0, L] x [0, T] K 2 : x G L°°(0, T; H 2 (0, L)), X G £ 2 (0, T; L 2 (0, £))}, 

endowed with the norm 

i 1/2 

(2-5) Hxllx : = esssu Pllx(-^)ll^ (0ii 

0<t<T 

which makes it a Banach space. It follows from the definition that 

(2.6) X c H 1 ^, T; L 2 (0, L)) with continuous embedding. 

Since every function x hi 1 (0, T; L 2 (0, L)) can be modified on a negligible subset of [0, T] 
so that i h-> xi't t) is strongly continuous from [0, T] into £ 2 (0, L) , we shall always refer to 
this modified function when we consider the properties of x(">*) f° r some i € [0,T]. With 
this convention we have 

(2.7) X c C°([0, T]; X 2 (0, L)) with continuous embedding. 
The following proposition shows the main properties of the space X . 
Proposition 2.2. Let x G X. Then for every t £ [0, T] we have x("i*) G H 2 (0,L) and 
(2-8) llx(^)ll ff2(0 ,L) ^ llxllx- 

Moreover, the function 1 1-> x("j *) is continuous with respect to the weak topology of H 2 (0, L). 
Finally, 

(2.9) xGG^M^aCL])), 
( 2 - 1() ) llxll c o([o,T];Ci([o,L])) < c WxWx > 

where the constant C is independent of x- 

Proof. To prove the first claim, let us fix t £ [0, T] and let N be a zero measure set up to 
which the essential supremum in ( I2.5D is actually a supremum. Consider a sequence t„ ^ N 
converging to t, so that \\x{;tn)\\ H ^ ,L) ^ IMIx • Since xMn) -> x(-^) in £ 2 ([0,£]) by 
02.7ft . we have that xG, *) G H 2 {Q, L) . Moreover, since the H 2 norm is lower-semicontinuous, 
we have also 

\\x(;t)\\ H z(a,L) <^v£\\x(-,tn)\\ H z(o,L) < llxllx. 

which proves 02.8D . Thanks to this inequality, the strong continuity of the function t i-> 
X (-,t) in L 2 ([0,i]) implies the weak continuity in H 2 (0,L). Then the compact embedding 
of H 2 (0,L) into C 1 ([0, L]) implies that the function t i-> x( - ,0 is continuous with respect 
to the strong topology of C^QO, L]) , which gives 02.9D . Finally, 02.1OI I follows from 02. 8D and 
from the continuous embedding of H 2 (0, L) into C 1 ([0, L]). □ 

Note that, by 02.9D . for every xe^we have 

(2.11) X,X' eC°([0,L]x[0,T]). 

We are interested only in functions X £ X such that s n> x( s j t) is the arc length parametriza- 
tion of a curve; this leads to the following definition 

(2.12) Ii:={x6X: \x'\ = 1 in [0, L] x [0, T] }. 
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3. Equations of motion 

In this section we derive the equations of motion for the swimmer. It is convenient to 
factorize the function x £ ^1 as the composition of a time dependent rigid motion r , which 
represents the change of location, with a function £ G X\ , which represents the change of 
shape. We write 

(3.1) x(s,t)=x(t)+R(t)£(s,t), 

where x(t) G M 2 is the translation vector and R(t) G R 2x2 is the rotation corresponding to 
the rigid motion r(i) . 

If we assume that / L £(s,i)ds = for every t G [O.T], then x(t) coincides with the 
barycenter of the curve which describes the swimmer at time t with respect to the 

absolute reference system, while the function £(-,t) will be regarded as the deformation 
seen by an observer moving with barycenter of the swimmer. 

Proposition 3.1. Let £ G X\ and let x: [0,T] M 2 and R: [0,T] -> IR 2x2 be functions such 
that R(t) is a rotation for every t e [0, T] . Then the following properties are equivalent: 

(i) the function x defined by 03.1D belongs to X\ ; 

(ii) the functions x and R belong to H r {Q,T). 

Proof. (i)=>(ii). For every t G [0, T] we define 

£(*) ~ \ J q £(*> *) ds and *(t) := | j£ xO, *) &s. 

Since {jeliC ff^O, T; i 2 (0, L)) we have that (je iJ^O, T) . By averaging CD} with 
respect to s we obtain 

(3.2) = x(t) + R(t)£(t). 
Subtracting this equation from d3.lt we obtain 

(3.3) X (M) -x(t) = R(t)(Z(s,t) -£(*)). 

Let us fix < G [0, T]. Since |£'(s, to) | = 1 for every s e [0,L], there exists s Q G [0, L] such 
that £(s ,io) — £(io) 7^ 0. By the continuity of £ there exist e > 0, an open neighborhood 
U of so G [0,L], and an open neighborhood V of to in [0,T] such that |x(s, t) — x(*)l = 
|£(s,t) — £(t)| ^ e for every s € U and every t G V , where the equality follows from 03.3D . 

Let X *{s,t) := ( X (s,t)-x(t))/ \x(s,t)-x(t)\ and £*(M) := (£(*,<)-£(*))/ |f(M) - • By 
03.31 1 we have 

X *{t,s)=R(t)C(t,s) 
for every s G U and every By elementary Linear Algebra we have 

( X *(s,t),e(s,t)) -( X *(s,t),JZ*(s,t)) ' 
(x*(s,t),Jti*(s,t)) { X *(s,t),Z*(s,t)) 

By (US and the functions x* and belong to H l {V; L 2 (U)) n C°(F; C^C/)) , so that 
the entries of the matrix in ( 13.41 ) belong to H 1 {V; L 2 {U)) . Since the matrix does not depend 
on s, we obtain R G H 1 ^). The conclusion R e H l (Q,T) follows now from a covering 
argument. 

Since x, R, and £ belong to H l {{), T) , we deduce from (US that x G iT^O, T) . 

(ii)=^(i). This implication follows easily from 02.41 1 and 03.lt . □ 

Remark 3.2. The purpose of the function £(-,t) is to describe the shape of the swimmer 
as a function of time. For each t we can choose the most convenient reference system. Of 
course, different choices are compensated by different rigid motions in 03.lt . 

In many cases it is convenient to describe the shape of the swimmer by means of the 
(oriented) curvature k(s, t) of the curve £(•, t) at s. This is because both in living organisms 
and in technological devices shape changes are usually obtained by controlling the mutual 
distance of several pairs of points. Prescribing the curvature can be interpreted as the 
infinitesimal version of this control, whose description is easier from the mathematical 
point of view. 

If x £ X\ and £ g X\ are linked by 03.1D . then clearly their curvatures are the same. 
Given £ G X\ , let #(s, t) be the oriented angle between the x\ -axis and the oriented tangent 
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to the curve £(•,£) at s. It is well known that «(s,t) = (£"(s,t), J£'(s,i)) — tf'(s,t), so we 
can easily get k from £ by differentiation and $ from ft by integration. In particular, if we 
assume £'(0, t) = e\ , we have ??(0, t) = 0, hence 



t?(s,i) 



Then the definition of #(s, i) gives that £'(s, i) 
we have 



k(ct, i) dcr. 

= (cos tf(s,i), sin tf(s,i)), so that, if £(0,f) 



£(M) 



(cos t?(<7, i), sin j?(cr, i)) dcr. 



This shows that the descriptions of the shape given by £(s, t) and k(s, t) are equivalent. 

By the change of reference 03.11 *. it is possible to rephrase the force and torque balance 
02.31 * and eventually obtain ordinary differential equations governing the time evolution of 
x and R . Those will be the equations of motion of the swimmer. We can write 



R(t) 



cos0(i) -sin0(i) 
sin0(i) cos0(i) 



where 9{t) is the angle of rotation. We assume that £ and x belong to X\ . Thanks to 
Proposition 13.11 we can differentiate 03. ID with respect to time. Plugging all the terms in 
0231 * and noticing that K x (s, t) = R(t)K s (s, t)R T (t) , we obtain 



(3.5) 



( m V 


' R{t) " 


{-[ 




1 





A(t) b(t) 
b T (t) c(t) 



R T (t) 





x(t) 
9(t) 



F sh (t) 
M sh (t) 



where M x (t) := M(t) - (F(t),Jx(t)) and H(t) 
matrix of 1161 . whose entries are given by 



A(t) b(t) 
b T (t) c(t) 



is the grand resistance 



(3.6a) 



A(t) := / K ( (s,t)ds, b(t) 



Kt(s,t)J£(s,t) ds, 



(3.6b) 



c(t) 



(J^(s,t),K € (s,t)J^(s,t)) ds. 



It is easy to see that the functions A, b, and c are ultimately determined by the shape 
function £ alone. The terms 



(3.7) 



F sh {t) :-- 



K € (s,t)i(s,t)ds, M sh (t) 



{JS(8,t),K 6 (8,t)i(8,t))d8, 



are the contributions to the force and torque due to the shape deformation of the swimmer, 
and they depend linearly on the time derivative £ . 

Enforcing the force and torque balance 02.31 * is equivalent to setting 03.51 * equal to zero 
and solving for x and , which eventually leads to the equations 



(3.8) 

where 
(3.9) 



x(t) 
0{t) 



R(t)v(t), 



v(t) := A(t)F sh (t) + b(t)M sh (t), u(t) := b T (t)F sh (t) + c(t)M sh (t), 

and Alt) , bit) , and clt) are the block elements of the inverse matrix 7^ _1 (t) . The structure 
of this system of ordinary differential equations is the same as that previously obtained in 
||3l|9j. The following result, analogous to [9, Theorem 6.4], holds 

Theorem 3.3. Let £ 6 X\ , x a e M 2 , and O G K. Then the equations of motion 03. 8D . with 
initial conditions x(0) = xq and 0(0) = 9 , have a unique absolutely continuous solution 1 1-> 
(x(t),6(t)) definedin [0,T] with values in R 2 xR. This solution actually belongs to H 1 (0,T). 
In other words, there exists a unique rigid motion t m> r(t), such that the deformation 
function defined by 03. ID belongs to X\, satisfies the equations of motion 02.3D . and the initial 
conditions x(0) = x and R(0) = Re , the rotation of angle O . 
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Proof. The result easily follows from the classical theory of ordinary differential equations, 
see, e.g., II15II . Indeed, the coefficients b T and c are continuous function of t, since they 
come from the inversion of the grand resistance matrix 1Z , whose entries are continuous in 
t . On the contrary, F sh and M sh are only L 2 functions of time. This is enough to integrate 
the second equation in 03.8D . By plugging the solution for 9 into the first equation and by 
an analogous argument on the coefficients A and b , also the equation for x has a unique 
solution with prescribed initial data. 

The last statement follows esily from Proposition [3J] □ 

Some notes on the matrix K and on the coefficients C T and C v are in order. First, we 
assume that 

(3.10) 0<C T <C v , 

secondly, we notice that the matrix K x (and therefore ) is symmetric and positive def- 
inite, and defines a scalar product in the space X\ . Indeed, by introducing the power 
expended during the motion 

V{x):=[ [ (-f(s,t), X (s,t))dsdt= [ f (K x (s,t)x(s,t), X (s,t)) dsdt 

(3.11) 



JO 
T r L 



[C T g(s,t) + C vX l(s,t)]dsdt, 



we find that 



Ct IIxIIl 2 (o,t ; l 2 (o,l)) ^ ^(x) ^ C v ||xlli,2(o,T;L 2 (o,L)) ■ 

Moreover, it follows from 02.2D and 02. IIP that the matrices and K x are continuous in 

(s,t). 

Finally, the strict inequality assumption C T < C„ cannot be weakened. Indeed, if we 
had C T = C v , then K x (s, t) would be a multiple of the identity matrix and therefore, from 
02.3al >. we would have 

= F(t) = -C T [ x(s,t)ds = 




JO ul - \ JO 

which is expressing that the barycenter does not move as time evolves. 

3.1. The shape function. We introduce now an important assumption on the shape func- 
tion £, called two disks condition, which rules out self-intersections of the swimmer. This 
hypothesis will be crucial in the proof of the existence of an optimal stroke. The idea un- 
derlying this condition is that two distinct points of the swimmer cannot become too close 
to each other during the motion. 

Definition 3.4. We say that £ e H 2 (0, L;R 2 ) with |£'(s)| = 1 for every s e [0,L] satisfies the 
two disks condition with radius p > if the following conditions are satisfied (see Fig. 0: 

(a) for every s e [0,L] there exist open disks B 1 (s), B 2 (s) of radius p > such that 
B 1 {s)nB 2 {s) = 0, £(s) eB 1 (s)nB 2 (s),andC((T) i B 1 (s)UB 2 (s) for every ere [0,L]; 

(b) there exist open half disks B~,B + of radius 2p centered at £(0) and £(L), respec- 
tively, with diameters normal to £'(0) and £'(£), respectively, such that £(<r) ^ 
B~ U B + for every a E [0, L] . 

Since £ is of class C 1 and |£'(s)| = 1, the disks considered in condition (a) are uniquely 
determined by £(s) and £'(s). Indeed, they are the disks with centers £(s) ± pJ£'{s) and 
radius p . In the sequel we will always assume that 

(3.12) B 1 (s) = B p ($(s) + pJl;'(s)), B 2 (s) = B p (£(s) - P J?{s)). 

The following proposition proves an important consequence of the two disks condition. 

Proposition 3.5. Let £ 6 H 2 (0,L;M. 2 ) with |£'(s)| = 1 for every s G [0,L]. Assume that £ 
satisfies the two disks condition for some radius p > 0. Then £ is infective on [0, L] . 
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Proof. Assume by contradiction that there exist si < s 2 such that £( s i) = £(s 2 ). It is easy to 
see that the two disks condition implies that £'(si) = ±£'(s 2 ). Assume that £'( s i) = £'(52) , 
the other case being analogous. Since these derivatives have norm 1, by changing the 
coordinate system we may assume that £'(si) = £'(s 2 ) = ei , the first vector of the canonical 
basis. We denote the coordinates of £ by £1 and £ 2 and we set a := £i(si) = £i(s 2 ) and 
P '■= £2(51) = £2(52)- By the Local Inversion Theorem, there exist e > 0, 5 > 0, and a 
C 1 function g : (a — e, a + e) — > R such that for every s e (s± - 5, si + (5) n [0, L] we have 
|£i (s) - a| < e and £ 2 (s) = Let 

U &{s), E '--= U ^OO- 

I s — s 1 1 <<5 |s — «i|<5 

By 03.121 1 it is easy to see that there exists an open rectangle R centered at (a, /3) such 
that E 1 n R = {(a,b) e R : b > g(a)} and E 2 n R = {(a, 6) e R : b < g(a)}. By condition 
(a) of Definition E3l£(s) G {(a, 6) E R : b = g(a)} for every s G [0,L] such that £(s,i) £ 
Therefore, s h> ((s) is locally an arc length parametrization of the graph of g . Since £(si) = 
£(s 2 ) = (a, 0) and £'(si) = £'(s2) = &\ , there exists 77 > such that £(si + s) = £(s 2 + s) for 
\s\ < 7], provided < sx + s < s 2 + s < L. This implies that £(cr— (s 2 — si)) = for every er 
in a neighborhood of s 2 in [0, L] such that a - (s 2 - si) G [0, L] . By taking the supremum a 
over c, we obtain that £(o- ) = for a Q := L- (s 2 — si) G [0, X) . If a > 0, we deduce also 
that £'(cr ) = • The same equality holds when a = 0, because in this case a = si = 
and s 2 = L, so that the equality follows from the assumption £'(si) = £'(s 2 ) • 

Let B + be the half disk considered in condition (b) of Definition 13.41 By the previous 
equalities, £(o"o) is the center of B + and £'(er ) points towards the interior. It follows that 
£(cr) G B + for some a > ct , and this contradicts condition (b). □ 



B 




FIGURE 1. Two disks condition: B 1 (s) and B 2 (s) are the open disks of 
radius p, B~ and B + are the open half disks of radius 2p. The dashed line 
represents the image of the map h defined in < I3.13I >. 



Given p > we introduce Cl, p ■= ([0, L]x{0}) + B p (0), the cigar-like set obtained by 
enlarging [0, L] x {0} . Define now a map h : [0, L] x (-p, p) — > R 2 by 

(3.13) h(8,v):=€(8)+yJ?(s). 

This map is extended to a continuous map /i: Cl, p — > R 2 by defining it as an isometry 
mapping Cl, p n {s < 0} into the half disk \B~ homothetic to the half disk B~ considered 
in condition (b) of Definition l3.41 with the same center and half the radius; the definition in 
Cl, p n {s > L} is similar and uses the half disk \B + (see Fig.Q}. The following proposition 
improves Proposition l3.5l and provides an equivalent formulation of Definition 13 .41 

Proposition 3.6. Let £ G iJ 2 (0,L;R 2 ) with |£'(s) = 1 for every s e [0,L], and let p > 0. 
The map h is injective if and only if £ satisfies the two disks condition with radius p. 

Proof. Let us assume that the two disks condition holds and let us consider two points 

(si,yi) ^ (s 2 , 2/2) in Gl, p ■ If < s\ = s 2 < L, then it must be \y± — y 2 | > 0, and therefore 
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|ft(si,j/i) - h(s 1: y 2 )\ = 1 2/1 - 2/2I > 0. If -p < si = s 2 < 0, then h(s 1 ,y 1 ) ^ h(s 1 ,y 2 ) since /i 
is an isometry on Cl. p n {s < 0} . The same conclusion holds if L < si = s 2 < L + p. 

Assume now that < si < s 2 < £• If 2/i = 2/2 = 0, then ft(si,2/i) = £(si) ^ £( s 2) = 
h(s 2 ,y 2 ), where the inequality follows from the fact that the curve is injective by 

Proposition 13.51 If 2/1 7^ = y 2 , then h(si,yi) belongs to one of the disks B 1 (si) ,B 2 (si) 
introduced in condition (a) of Definition 13.41 It follows that h(s 2 ,y 2 ) = £(s 2 ) ^ B 1 (si) U 
B 2 (si) , hence h(si,yi) ^h(s 2 ,y 2 ). The same conclusion holds if 2/1 = 7^ y 2 ■ 

Let us consider now the case < si < s 2 ^ L and 2/1 7^ 7^ 2/2 • Define := {£(s,) + 
yJ£'(si) : < (sign 2/1)2/ < p} for i = 1,2. Let us prove that 

(3.14) St n S 2 = 0. 

Let p,; := £(si) and let A be the open disk with center q := £(s») + sign(i/j)p JC( s i) an d 
radius p. Note that and pj are the endpoints of S 1 , and that pi 7^ p 2 since the map 
s \-> £(s) is injective by Proposition 13.51 If ci = c 2 , then 5i n 5 2 = because 5i,52 are 
radii of the same circle with different endpoints. 

We consider now the case a 7^ c 2 . Since Di is one of the disks B 1 (s±), B 2 (si) , by 
condition (a) of Definition 13.41 we have that pi = £(si) 6 dDi and p 2 = £(s 2 ) ^ £>i ■ 
Similarly, we prove that p 2 e dD 2 and pi D 2 . Therefore, pi G dD{\D 2 and p 2 G dD 2 \D 1 . 

Assume by contradiction that Si and S 2 meet at some point p, which must belong to 
Di n D 2 . Let z be the intersection between dD 2 and the half-line stemming from c 2 and 
containing ci ; under our assumptions, we have z G Di . Since p e Di n D 2 and pi G 
dDi \ D 2 , there exists a unique point g G dD 2 n Di on the segment joining p and pi . Now, 
the half-line through p stemming from c 2 meets dD 2 on the smallest arc T with endpoints 
q and z. Since q, z G Di and the disks have the same radius, we have rcDi (see Fig.©. 
The previous argument shows that p 2 G T, which contradicts the condition p 2 G <9Z? 2 \ Di . 
This concludes the proof of the equality Si n 5 2 = in the case < si < s 2 ^ L and 
2/i 7^ 7^2/2 , and implies that h(si,yi) ^ h(s 2 ,y 2 ). 

Ji 

^ V \ 





FIGURE 2. Injectivity for < si < s 2 ^ L and 2/1 ^ 7^ j/ 2 in the case 
ci 7^ c 2 : Two possible situations contradicting 03.14D . 

We consider now the case < si < L < s 2 < L + p. Assume by contradiction that 
h(si,yi) = h(s 2 ,y 2 ) =: p. Observe that p G \B + . Denote pi := £(si) and S := {£(£) + 
yJ£'(L) : —p < y < p) . By 03.14D for si and L the segment with endpoints p,pi does not 
intersect S. On the other hand |p-pi| = \yi\ < p. By elementary geometric arguments 
we find that the set Q of points which can be connected to a point of \B+ by a segment 
disjoint from S and of length less than p is contained in the union B + U B 1 (L) L) B 2 (L) (See 
Fig. [3). Therefore pi = £(si) G Q and this violates either condition (a) or condition (b) in 
Definitionl3"74"l 

In the case si = L < s 2 < L + p, we have h(si,yi) G d(\B + ), while h(s 2 ,y 2 ) G \B + , 
so that h(si,yi) ^ h(s 2l y 2 ). The cases — p < si < < s 2 < L and — p < si < = s 2 are 
analogous. 

The last case to consider is when —p < si < and L < s 2 < L + p. Assume, by 
contradiction, that h(si,yi) = h(s 2 ,y 2 ). Since the radius of curvature of £ is always less 
than p, one can prove (see Lemma 13771 below and Fig. [4]) that 

(3.15) D-cB-U |J (B 1 (s)U{^(s)}UB 2 (s)), 

0<s<f P 
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Figure 3. Injectivity for < si < L < s 2 < L + p. The shaded region 
represents the set Q ; the thick segment is the set S . 

where D~ := {x G M 2 : d(x, \B~)< p} and B 1 (s), B 2 (s) are the open disks defined in 03.12D . 




FIGURE 4. The region D~ of 03.15b . 



Since h(si,yi) G \B , h(si,y\) = h(s 2 , 1/2), and |/i(s 2 ,y 2 ) — £(L)\ < p, we deduce that 
£(£) e D-. By dBTTSl either £(L) G B~ or there exists s G [0, f] such that £(X) G S x (s) U 
{£(s)} U B 2 (s). This contradicts Definition 13.41 or Proposition [33] and concludes the proof 
in this case. 

Let us now assume that h is injective, and consider s G [0, i] . Let B 1 (s), B 2 (s) the open 
disks defined in (EH). It is clear that B 1 {s)nB 2 (s) = and B 1 (s)nB 2 (s) = {£(s)}. Denote 
by N(s) := {S,(s) + yj£'(s), y G K} the normal line through the point £(s), and let £(s) := 
£(s)+-^h\Jt;'(s) be the evolute of £, i.e., the curve that contains the centers of the osculating 
circles to £. Since £ is a curve of class H 2 , the curvature is well defined almost everywhere; 
let s\ e (0,i) be a point at which the curvature k(si) is defined, and let s 2 G (0, L) be 
another point. By the injectivity of h, the normal lines N(si) and -/V(s 2 ) cannot meet at a 
distance less than p from the curve £, and their intersection tend to £(si) as s 2 — > si (see 
IfTTl Ex.7, page 23]). This implies that k(si) < 1/p. We conclude that k(s) < 1/p for a.e. 
s e [0, L] . Let I p (s) := {a e [0, L] : \cr — s\ < p} n [0, L] be a neighborhood of s of radius p in 
[0,L], By standard results in Differential Geometry, £(cr) B 1 (s)L)B 2 (s) for every a G I P {s). 
If now a G [0, L] \ I p (s) , then again £(er) ^ ^(s) u B 2 (s) . Indeed, assume by contradiction 
that there exists a € [0,L] \ I p (s) such that £(cr) G B 1 (s). Then there exists r G (0,p) such 
that^(a) € B r (Z(s)+rJ£'(s)). Then, let Vl := min CT g [0>L] {|C(s) + rJ£'(s) -£(ct)|} < r, and let 
si G [0, i] \ I p (s) be the point where the minimum is attained. If si G (0, L) , then 

either + rJ£'(s) = £( Sl ) + yi J^(si) or £(s) + r J£'(s) = £( Sl ) - yi J£'(si). 

This violates the injectivity of h. The cases s\ = and s\ = L lead to a similar contradic- 
tion, taking into account the definition of h near the endpoints of the segment. This proves 
that condition (a) in Definition ^. 4l holds. 

To prove that condition (b) holds, we assume, by contradiction, that there exists s G 
[0,L] such that £( s o) G B~ . We first observe that, for ,s g [0,717?], the point £(s) lies on 
the opposite side of B~ with respect to its diameter {£(0) + yJ£'(0) : —p < y < p] , so that 
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£(s) B for s G [0, irp] . Then s belongs be the closed set £ of points s e [7173, £] such that 
£(s) lies in the closure of B~ . Let s\ be the minimum point of 

min|£( S )-£(0)|. 

Since £(s) ^ S^O) U £? 2 (0) for every s <E [0,L], the point £(si) does not belong to the 
diameter {£(0) + yJ£'(0) : —p <y < p}. Therefore £(si) e B~ and si belongs to the interior 
of S. This implies that £(si) — £(0) is orthogonal to £'(si), hence £( s i) — £(0) = yJ£'(si) for 
some y G (-2p, 2p) . Then ft.(si, — |) = £(si) - | J^'(si) G |-B~ . This violates the injectivity 
of h and concludes the proof of the proposition. □ 

Lemma 3.7. Let p > and let £ € H 2 (0, f ,o;K 2 ) ivii/i |£'(s)| = 1 f° r every s e [0, § p] and 
|£"(s)| < l//o /"or a.e. s e [0, f p]. TTie/i < f37T5l > holds. 

Proof. Let £i(s) and £ 2 (.s) be the coordinates of £(s) and let i9(.s) be the oriented angle 
between ei := (1,0) and £'(s). It is well known that |$'(s)| = |£"(s)| < I// - Assume, for 
simplicity, that £(0) = and £'(0) = ei . By standard results in Differential Geometry the 
curve £ does not intersect the open disks B p (±/9e 2 ), where e 2 := (0, 1). By integration we 
obtain 

£i( s ) = / cosi9(cr) der . 
Jo 

Let us prove that 

(3.16) psm(s/p)<^(s) 

for every ,s e [0, |]. Since the inequality is true for s = 0, it is enough to show that 
cos(s/p) < ££(s) = cost?(s) = cos|$(s)| for every s e [0, §]. Since £'(s) = (cos$(s),sin$(s)), 
we have |$'(s)| < 1/p, hence |$(s)| < s/p, so that cos(s/p) < cos |i?(s) by the monotonicity 
of cos in [0, §] . This concludes the proof of ( 13.16t . 

Let S be the segment with endpoints (p,±p), which belong to the circles dB p {±pe 2 ) ■ 
Inequality ( 13.161 * implies that the curve £ intersects the segment S . Since the bound on 
the curvature implies that £ cannot have a vertical tangent, except when £ is contained 
in dB p (±pe 2 ) , the intersection point is unique. Let s e [0,p-|] be the value of the arc 
length parameter of this intersection point and let N(s ) :~ {£(so) + yJC( s o),y £ ^} be the 
corresponding normal line to £ . 

If £(.so) = (p,±p), then the bound on the curvature implies that £ is contained in 
dB p (±pe 2 ) and the statement of the lemma is easily checked. So we may assume that 
£( s o) (/°! ±p) ■ We may also assume that C 2 ( s o) > 0. If not, we just reverse the orientation 
of the x 2 -axis. 

Let B 1 (s ) be the tangent disk to £ at £(so) denned by B 1 (s ) := B p (£(sq) + pJ£'(s )). 
We claim that 

(3.17) (p,p) €B\s ). 

If £ 2 (s ) = 0, then (p,p) e N(s ) and its distance from £(s ) is less than lcngth(S) = 2p, 
which implies 03.17D . If £ 2 (s ) =/= 0, we argue by contradiction. Assume that 03.17D is not sat- 
isfied. Then dB 1 (s ) intersects S in £(s ) and in another point £ between £(s ) and (p,p). 
Therefore the center c\ of the disk B 1 (s ) is the vertex of an isosceles triangle with basis 
contained in S and equal sides of length p . Elementary geometric arguments show that 
this vertex must belong to the astroid obtained by removing the four disks B p ((±p, ±p)) 
from the square (— p, p) x (— p, p) (see Fig. 0. Therefore the distance from the origin of the 
center of B 1 (sq) is less than p. This implies that e B 1 (s ). On the other hand, since £ 
is tangent to this disk at £(so), the bound on the curvature implies that £(s) ^ B 1 (s ) for 
every s e [0, §] . This contradicts the assumption £(0) = and concludes the proof of 03.17D . 

Let C be the curvilinear triangle obtained by removing the disks B p (±pe 2 ) from the 
rectangle (0,p) x (— p, p). Let C+ and C~ be the parts of C weakly above and strictly 
below N(s ) and let p be the intersection of N(s ) and <9£? p (pe 2 ) contained in the closure 
of C . Since the distance between p and £(so) is less than 2p, we deduce that p e B 1 (s ). 
Since (p, p) e B 1 (s ) and £(s ) € dB 1 (s ) , we obtain that C+ is contained in i? 1 (s ) . 

Let us prove that 

C-c |J (S 1 ( S )u{e( S )}U J B 2 ( S )) 

0<s<s 
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FIGURE 5. The absurd situation in the proof by contradiction. The circle is 
B 1 (s ) , while the curve is not shown. 



Let us fix x e C and let s\ be a minimizer of 

min \x-£(s)\. 

0<s<so 

Since (x,£'(0)) > and (x - £(so), C'( s o)) < 0, we have si e (0, s ), hence the orthogonality 
condition 

(x-e(si),r(si)) = o. 

Since \x— £(si)| < \x\ < 2p, the point x belongs to B 1 (si)U{t;(si)}UB 2 (si) , which concludes 
the proof. □ 

We now prove a result stating that a bound on the angle ■& formed by the tangent with 
the first axis implies the non self-intersection of the swimmer. 

Lemma 3.8. Let ?9 e C^QO,!,]), let k := max{|i9'(s)| : s e [0,L]}, and let 

Assume |$(s)| < 7r/4 for every s e [0,L]. Then, £ satisfies the two disks condition with 
radius p, for every < p ^ 1/k . 

Proof. Notice, in the first place, that |£'(s)| = 1 for every s e [0, L] , so that £ is a regular 
curve parametrized by arc length. The condition \"&{s)\ < 7r/4 for all s e [0, L] implies that 
£ is a graph with respect to the xi-axis. Given p, with < p ^ 1/ko, we define the open 
disks B 1 (s) := B p (£(s) + pJ£'(s)) and B 2 {s) := B p (£(s) - pJ£'(s)), as in ( 13321 Since £ is a 
graph and its curvature is bounded by k , the disks i? 1 (s) and B 2 (s) satisfy condition (a) 
in Definition 13.41 Finally, the construction of B~ and B + as in condition (b) of Definition 
13.41 is straightforward. □ 

The preceding lemma will be useful in Section [4] to check that the deformations we con- 
struct to prove the controllability of the swimmer are admissible. 

4. Controllability 

In this Section we show that the swimmer is controllable, i.e., it is possible to prescribe a 
self-propelled motion that takes it from a given initial state Xin to a given final state xfin • 
More precisely, we prove the following theorem. 

Theorem 4.1. Let p > and let Xm, Xfin £ H 2 {Q,L;R 2 ), with |x( n (s)| = |Xfi n 0)l = 1 for 
every s £ [0, L] . Assume that Xin and Xfm satisfy the two disks condition with radius p 
(see Defi.nition \3.4h Then, there exists x 6 Xi , satisfying the force and torque balance 02. 3D . 
such that xi s > 0) = Xin(s) and x( s , T) = Xfm(s) for every s e [0, L], and such that for every 
t E [0, T] the curve x(-, t) satisfies the two disks condition with radius p. 

Proof. To construct x> we divide the interval [0,T] into three intervals [0, |T], [|T, |T], 
and [|T, T]. In the first interval we straighten xin , i-e., we construct x> satisfying the 
force and torque balance 02.3D and the two disks condition with radius p on [0, |T] , such 
that x(sj 0) = Xin(s) and x( s > = E in (s) for every s e [0, L] , where S in is the arc length 
parametrization of a suitable segment of length L , depending on Xin • 



ONE-DIMENSIONAL SWIMMERS IN VISCOUS FLUIDS 



13 



The same construction, with time reversed, shows that there exists a segment Eg^ , de- 
pending on Xfin , that can be transferred onto xtm > Le., there exists x satisfying the force 
and torque balance 02.3D and the the two disks condition with radius p on [|T, T] , such that 
X(s, |T) = E fln (s) and x(s, T) = Xfin(s) for every s e [0, L] . 

Since, in general, Efm does not coincide with E in , we use the interval [|T, |T] to transfer 
E in onto S fin . 

We now describe the construction of x on [0, |T] . First of all, it is possible to find S, £ Xi 
such that £(s,0) = Xm(s) for every s e [0,L] and s >-> £(s, ^T) is affine on [0, L], It is also 
possible to obtain that satisfies the two disks condition with radius p > for every 

ie[0,±T]. 

The last requirement can be fulfilled in the following way. If at one end of the swimmer 
there is enough room, we pull it along the tangent and unwind it from its original shape 
obtaining a straight configuration, as illustrated in Fig. [6] 




Figure 6. Straightening the swimmer I. The dashed line is the straight- 
ened configuration. 

If this is not the case, then we operate as in Fig.[7j the unwinding is achieved by pinching 
a point with maximal x\ -coordinate and pulling it to the right respecting the two disks 
condition. 




Figure 7. Straightening the swimmer II. The dashed lines represent the 
intermediate phases of the stretching procedure. 

We now compose £ with a time dependent rigid motion and define x on [0, L] x [0, |T] 
by 03.lt . Clearly the curve x(-, t) continues to satisfy the two disks condition with radius p 
for every t e [0, ^T] . Moreover the function Si n (s) := x( s i is affine and E[ n (s)| = 1 for 
every s e [0, L] . The vector x and the rotation R are chosen so that the equation of motion 
( 13.81 1 is satisfied (this is possible thanks to Theorem l3.3D . so that x satisfies the force and 
torque balance 02.3D in [0, \T] . 

Note that, while the affine map £(•, \T) can be chosen freely, the corresponding map 
S in depends on the superimposed rigid motion, which, in turn, depends on the data of 
the problem. Therefore, in this construction the location of the segment S in cannot be 
prescribed. 

On T] the function x is defined in a similar way. To transfer S; n into E nn in the 
time interval |T] we show that, for a straight swimmer, it is possible to produce self- 
propelled motions achieving any prescribed translation along its axis and any prescribed 
rotation about its barycenter. 

To summarize, the whole control process is organized as in Fig. [8) 

, . straightening . . rotation, translation, rotation , , straightening ~ 1 , , 

Xin(0 > Ein(-) > Efm(-) > Xfin(-)- 
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''Si 



in 



Figure 8. Sketch of the control process. 



4.1. Translation. In this subsection we describe how to translate a straight swimmer in 
a given time interval: since the problem is rate independent (see Remark 12. It . it is not 
restrictive to work in the interval [0, 1] . The motion of the swimmer is obtained through 
the translation along the swimmer itself of a localized bump. In order to get a rectilinear 
motion of the endpoints of the segment, we have to assume that the bump satisfies some 
symmetry properties. 

We can distinguish three phases of this motion, which now takes place in the conven- 
tional time span [0, 1] : the transient of formation of the bump [0, to] , with < t < 1/2, the 
translation of the bump [to, I— to] , and the transient of destruction of the bump [I— to, 1] - We 
can assume that the swimmer lies initially on the x\ -axis and that the initial parametriza- 
tion is Si n (s) = sei . As before, the motion will be first described through a function ^ e Xi 
satisfying the two disks condition with the prescribed radius p . Then, £ will be composed 
with a time dependent rigid motion in order to obtain \ satisfying also the force and torque 
balance. 

For simplicity, we assume £(0, t) = for every t £ [0, 1] . The function £(s, t) will be better 
described by means of the angle between its tangent line and the positive xi-axis. 

This leads to the formula 



To choose §{s,t), we fix a smooth function $q : K — > (— tt/4, 7t/4), with support spti? = 
[—1,1], a constant I e (0, L/2), and a decreasing affine function (: E — > K such that 

C([*o,l-*o]) c [I, L - I] ; for every (s,t) e [0,L]x[0,l] we set 



Thanks to Lemma [3~8l it is possible to choose the curve £(•,£) satisfies the two disks condi- 
tion with radius p for every t e [0, 1] . 

Note that for every t e [t , 1 - t ] we have spt «?(-, t) = [((t) - I, ((t) + 1} . Since §{s, t) = 
for s e [0,C(t) - I], we obtain from d4TTl > that ^(s,i) = se x , for s £ [0, C(<) - Similarly, 
since #(s, t) = for every s 6 [C(t) + Z, L], 04. ID implies that in this interval £{-,t) is the 
arc-length parametrization of a segment parallel to the x\ -axis. Therefore, the curve £(•, f) 
is the union of two segments and a connecting bump, corresponding to the restriction of the 
curve £(■, t) to the interval [((t) — I, + I] . Since ( is decreasing, the bump is moving to 
the left. 

Moreover, by a change of variables we have 



(4.1) 




(4.2) 





Together with the previous remarks, this implies that 



(4.3) 
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For every t 6 [0, to] , the support of $(•, t) is independent of t and coincides with [((t ) - 
h((to) + I], while the size of $(s,t) depends on t and tends to zero as t — > 0. The cor- 
responding function £(,s, t) given by 04.11 * can be interpreted as a transition from the seg- 
ment parametrized by sei , s G [((to) — ^C(^o) + /], to the bump parametrized by £(s, to), 

se [C(*o)-*,C(*o)+*]- 

Some symmetry assumptions will be made on the angle function t? which generates 
the bumps. We recall that a function u: [a, b] — > K is said to be even (resp. odd) in [a, b] if 
u(x) = u(a + b — x) (resp. u(x) = —u(a + b—x)) for every x e [a, b] . The main consequence of 
this symmetry properties is summarized in the following lemma, whose elementary proof 
is omitted. 

Lemma 4.2. Let u: [a,b] — > K be an even (resp. odd) function and let c := (a + b)/2 be the 
middle point of [a 7 b]. Then, the integral function U(x) := u(s)ds is odd (resp. even) in 

[a, b] . 

Given a e K , we describe a self-propelled motion that transfers the segment Ej n (s) = sei 
into the segment Sfi„(s) = (a + s)e\ , s e [0,L]. It is not restrictive to assume a > 0. 
We make the following symmetry assumption on the angle function t? : 

(4.4a) O is odd in [-1,1]; 

(4.4b) o is even in [-1, 0] and in [0,1]; 

(4.4c) z?o is odd in [-1,-1/2], [-1/2,0], [0, 1/2], and in [1/2,1]. 

Notice that 04.4aD yields /_ x sin$o(s) ds = 0, therefore, 04. ID and 04.2H imply that for 
every t e [to, 1 - *o] the curve s h-> t;(s, t), s e [C(*o) +l,L], parametrizes a segment lying on 
the xi-axis. 

Figure [9] shows an example of a bump £(•, t) , for t e (t , 1 — to) , whose angle function $ 
enjoys the properties listed above; notice that the parity of the vertical component of £ is 
reversed with respect to that of . 




FIGURE 9. Image of the function s h-> for t £ {t ,l - t Q ) for the 

translational motion. 



We now choose the affine function ((t) so that C(*o) = L — l and C(l —to) = l. This gives 
the expression 

C(t)=L-l-c(t-t ), 

with 

(4.5) c:= (L — 2l)/(l — 2t Q ). 

We now compose ^ with a time dependent rigid motion and define \ on [0, L] x [0, 1] by 
03.H . The vector x(t) and the rotation R(t) are chosen so that the force and torque balance 
( 12.31 1 is satisfied by x ■ 

We make the Ansatz that R(t) = I and x 2 (t) = for all t e [0, 1] . We will show that it is 
possible to determine xi(t) in such a way that x satisfies 02.3I I. 

With these assumptions, the expression for x reads x(s,t) = xi(t)ei + £(s,t). Since 
^(s, 0) = £(s, 1) = sei , the initial and final conditions for x mentioned after Lemma F4.2l are 
equivalent to 

(4.6) xi(0)=0, x 1 (l) = a. 

From the formula for x we get x'(s,t) = ^'(s,t) and x(s,i) = xi(t)ei + £(s,t). It follows that 
the matrix K x defined in 02. 2D satisfies 

(4.7) K x (s,t) = Kt(s,t) = C T t'(s,i)®Z'(s,t)+C u (Jl;'(s,t)) ® (J$'(s,t)). 
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The linear densities of force and moment (see 02.1D ) are then given by 

-f{s,t) = K x (s,t)x{s,t)=x 1 (t)K e (s,t)e 1 +K c (s,t)£(s,t) 
-m(s,t) = x 1 (t)(K i (s,t)ei, J£(s,t)) + ±i {t)x\ (t) (K^(s, t)e x , e 2 ) 
+ (K ( (s, t)£(s, t), J£(s, t)) + Xl (t)(K^s, t)£(s, t), e 2 ) 
By plugging this information in 02. 3D . we get 

(4.8) ±i(t) Kt(s,t)e 1 ds = - K s (s,t)i(s,t) ds. 

Jo Jo 

&i(t) J (K i {s,t)e 1 ,J£,(s,t))ds + x 1 (t)x 1 (t) j (K i (s,t)e 1 ,e 2 ) ds = 

(4.9) 



(Kt(s,t)ZM,Jfa,t))te-xi(t) / (K^s,t)^(s,t),e 2 ) ds 

Jo 

To solve simultaneously these equations for the unknown xi(t), we will show that the 
second components of the integrals in 04. 8D are zero, that the first component of the integral 
in the left-hand side in 04.8H is non zero, and that all integrals in 04.91 * are zero. 

4.1.1. Stationary regime. We now solve 04.8H and 04.9D for t e [to, 1 — to]. We first deal 
with the left-hand side of 04.8D . which can be written as 

f L Mt)-l K(t)+i pL 

(4.10) / K i (s,t)ds = / K i (s,t)ds + / K i (s,t)ds + / K s (s,t)ds. 

Jo Jo J((t)-l Jc(t)+i 

Since £'(s,t) = e t and J£'(s,t) = e 2 for s £ [((t) - l,((t) + I], by j4T7l > the off-diagonal 
elements of the first and third integrals in 04.1OD vanish. To show that the same is true for 
the second integral, we write it as 

sin 2 do — cos $o sin do 



IC T 



cos i?o cos $0 sin 
cos do sin do sin 2 do 



ds + IC V 



— cos do sin do cos do 



ds. 



Thanks to 04.4al > the function s t-> sin$ (s) cos$o(s) is odd in the interval [—1,1], therefore, 
all the off-diagonal elements in the matrices above integrate to zero. Similarly, we can 
prove that the diagonal elements of the integrals in O4.10l > are strictly positive. 

Let us show now that the second component of the integral in the right-hand side of 04.8D 
vanishes. First, notice that for s [((t) - l,((t) +!], we have £(s,t) — by 04.3H . so that 
the integral we have to study reduces to an integral on the interval between ((t) — I and 
C(t) + 1. Recalling that * e [to, 1 - to], by gj} and J4~2l ). for every s e [((t) - I, ((t) + 1} we 
have 



Hs,t) = j f ( -^ffif'f ) d> W(a,t))d* = c ( C ° S *ffr\- 1 ) 
lj«t)-l\ cosd {ip(<T,t)) J y smd (^{s,t)) J 

i>(s,t) 



>C(t)-i 

where 

s - C(t) 
I 

Here we have used the fact that the integrand above is a total derivative. Therefore, by 
changing variables, we obtain 

C(t)+i 

Kz(s,t)£{s,t) ds 

C(t)-; 

= dc r t ( cos tf J ) (( cos t°!1 ) , ( cosM q s ] ; 1 ))ds+ 

_ ± \ smdo{s) J \\ smdo(s) ) \ sini? (s) J/ 
1 / -smtfo(s) \/( -smdo(s) \ ( cosdo(s) - 1 ' ' ^ 



cos"i?o(s) /\\ cos ^o( s ) J\ sin do (s) 
To prove our claim we will show that 



(4.12) J cos d (s) sin d (s)ds = J sintf (s)ds = 
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By 04.4aD . it is easy to notice that the integrands above are odd functions in [—1,1], there- 
fore, 04.121 1 holds. This concludes the proof that the second component of the right-hand 
side of 04.8D is zero. Therefore, for every t e [to, 1 - to] we can solve ( 14.81 1 for x±(t) uniquely. 

We must now prove that also 04.91 1 is satisfied for every t e [to, 1 — to] . This equation 
reduces to 



(4.13) 



±i(f) / (i^(M)ei, J£{s,t))ds 



(K ( (s,t)i(s,t),JZ(s,t)}ds, 



since the off-diagonal elements of the integral in the left-hand side in 04.81 * are zero, as 
well as the second component of the integral in the right-hand side of the same equation. 
To prove that 04.13D is satisfied, we show that both integrals in 04.13D are zero for every 

te [t ,l-to]. 

To prove that the first integral is zero, by 04.7D it is enough to show that 



<£'(*> *), ci> (£'(*, t), J£(a, t)> ds = / (J£'(a, t), ei) (J?(s, t), J£(s, t)) ds = 0. 



(4.14) 



Since £'(s, i) = e% and J£(s, t) is proportional to e 2 , for every s £ [((t) - l,C(t) + l\> these 
integrals reduce to integrals on {((t) — l,C(t) + We can write the integrals in 04.141 ) 
explicitly using 04. ID and the second line in 04. 2D and we can change the variables as above. 
Recalling that the off-diagonal elements of 04.1OD are zero, in order to prove 04.14D we have 
to show that 



(4.15) 



cos do(s) sin do(s) 



cos $o(c) da- 



ds = 0, 



(4.16) 



(4.17) 



sin 2 i?o(s) 



-l 
l 



cos 2 #o(s) 



sin $o(a) da 



sin$o( cr ) da 



ds = 0, 



ds = 0. 



Since cos$ (s) is even and sin$ (s) is odd in [—1,1] by 04.4aD . we have f_ x cosz? (s) sini? (s) ds 
0. Therefore, to prove 04.15D it is enough to show that 



(4.18) 



cos i?o( s ) sin $o(s) 



cos $o(a) da- 



ds = 0. 



By Lemma 04.2D . the function s t-s> f Q s cos $o(a) da is odd in [—1,1], and consequently the 
integrand in 04.18D is even. Hence 04.18D reduces to 



(4.19) 



cos $o(s) sin do(s) 



cos $o(c) da 



ds = 0. 



Since cosz? (s) and sini9 (s) are even in [0, 1] by 04.4bD . we have 



1/2 

cos"#o(s) sini9 (s) ds = 2 / cos-do(s) sini5 (s) ds = 0, 



(4.20) 



where the last equality follows from the fact that cosi? (s) is even and sin$ (s) is odd in 
[0, 1/2] by 04.4cD . This equality implies that 04.19D reduces to 



(4.21) 



/' 

Jo 



cos do(s) sin flo(s) 



1/2 



cos $o(a) da 



ds = 0. 



Since cos-# (s) and sini9 (s) are even in [0, 1] by 04.4bD and s m> f*, 2 cosi9 (a) da is odd in 
[0, 1] by Lemma [4721 we obtain 04.21D . This concludes the proof of 04.15D . 
We now prove 04.16D . Since sin$ (a) is even in [— 1, 0] by 04.4bD . we have 



(4.22) 



sin 'do (a) da 



1/2 



sin^o(a) da = 0, 
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where the last equality follows from the fact sm<d (a) is odd in [-1/2,0] by 04.4cl >. There- 
fore, to prove 04.161 * it is enough to show that 

r-l 



(4.23) 



sin i?n(s) 



sin^o((7) dcr 



ds = 0. 



Since siin9 ( s ) is odd in [—1, 1] by 04.4aD and s m> f Q s sin $ (cr) dcr is even in [—1,1] by Lemma 
14.21 equality 04.231 * reduces to 

r-l 



(4.24) 



sin i?n(s) 



sin $n(cr) dcr 



ds = 0. 



1/2 

Since sm$ (c) is odd in [0,1/2] by 04.4cl * we have J Q sin -do(a) da = 0, hence 04.24D is 
equivalent to 



(4.25) 



/ sin 2 # n (s) / sin ^o(o-) do- 
Jo A/2 



ds = 0. 



Since sin^s) is even in [0, 1] by d4.4bl ) and the function s n> f*. 2 sin #o(ct) dcr is odd in [0, 1] 
by Lemma 14.21 equality 04.251 * holds. This concludes the proof of 04.161 *. Equality 04.17D is 
proved in the same way, thus the proof of 04.14D is complete, so that the left-hand side of 
04.131 * is zero. 

Similar arguments show that also the right-hand side of 04.131 * is zero. 
Therefore we have proved the existence of a function xi(t) such that x(s,f,) = xi(t)ei + 
£(s, t) satisfies the force and torque balance 02. 3D for every t e [to, 1 - t ] • 

4.1.2. Transient regimes. We now solve 04.8D and 04.9D for t e [0, to] , recalling that ((t a ) = 
L — I . We first deal with the left-hand side of 04. 8D . Arguing as in the stationary regime, we 
can prove that the off-diagonal elements of the matrix 

fL 

K^(s,t) ds 



are equal to zero. The only difference in the proof is that now, since t e [0, to] , 

cos(^oWs))) 



where 



sin(£tf (V(s))) 
s-L + l 



I 



We now prove that the second component of the integral in the right-hand side of 04.81 * is 
zero. As in the stationary case, the only contribution is on the interval [L — 21, L] , where 
£(s, t) is now given by 



-sin(±tf (V(cr))) 



h Jl-21 V COS^oWKO)) 

Therefore, after changing variables, we obtain 
(4.26) 



MM*))**- 



Ke(a,t)£(8,t) ds 



l 2 C T 



l 2 C v 
to 



cos(£t?oOO) 

M-&M*)) 

cos(-t-Ms)) 



cosies)) 
-sin(£0 o («)) 



i?oO)dcr)ds 



-sin(^ (cr)) \ \ 
i t o i w #o lender ) ds. 



To prove our claim we will show that 



(4.27) 
(4.28) 



cos(£0oOO)sin(£0oOO) 



sin 2 (^ (s)) 



sin(£0 o (<r))0o(<r) dcr 
cos(f ^o(o-))^o(ff)dcr ( 



ds = 0, 
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(4.29) 



cos 2 (^ (s)) 



c M (^oW)tfoWda 



ds = 0. 



Since cos(^-#o(s)) sin(^-i? (s)) is odd in [— 1, 1] by 04.4aD . equation 04.271 * is equivalent to 



(4.30) 



cos(£0„(s))sin(£0oOO) 



sin(£tf (<T))tf (a)d<7 



ds = 0. 



By 04.4al * and by Lemma 14.21 the function s h> J q s sin(^-i9o(c))^o(c) da is odd in [—1,1], 
therefore 04.301 * is equivalent to 



cos(£tf (s))sin(£tf ( s )) 



sm(±Mv))Mv)da 



ds = 0, 



since its integrand is even in [—1,1]. By 04.201 *. this is equivalent to 

»i 



(4.31) 



cos(£0oOO)sin(£0oOO) 



1/2 



ds = 0. 



By d4.4bl >. cos(^z9 (s)) , sin(^-$ (s)), and sin(^#o(s))#o(s) are even in [0, 1] , hence the func- 
tion s i y f*, 2 sm(-^-d (a))-do(a ) da is odd in [0,1] by Lemma [4721 This implies that 04.31D 
holds, since its integrand is odd in [0,1]. This concludes the proof of 04.27I *. 

To prove 04.28D we notice that in the function s h> cos(^#o(s))#o(s) is even in [—1,0] by 
04.4bl *. Hence, 



(4.32) 



cos{{- o d {a))M°) da = 2 



-1/2 



cos(£0o(<t))0o(<O da = 0, 



where the last equality follows from the fact that the function s h-» cos(^-$ (s))^o(s) is odd 
in [-1, -1/2] by 04.4cl *. This implies that 04.281 * is equivalent to 



(4.33) 



sin 2 (-^ (s)) 



cos{j-J {a))M<y)da 



ds = 0. 



By 04.4al *. the function cos(^-$o(c))$o(c) is odd in [—1, 1] , and therefore the function s n> 



/ cos(^-i?o(o'))i?o(o') dcr is even in [-1, 1] by Lemma I4TZ1 Since also sin (j-$o{s)) is even in 
[-1,1] by 04.4al *. 04.33D is equivalent to 



sm 2 (±Ms)) 



(i 



cos(£0 o (<r))l?o(<r)dff 



ds = 0, 



which, by 04.32D . is equivalent to 

/ 



(4.34) 



sin 2 (^ (s)) 



1/2 



ds = 0. 



By 04.4bl *. the function cos(^i? (o'))^o(c r ) is even in [0, 1], and therefore the function s n- 



^0 

Jy 2 cos(j^'3a(a)){)o(a) da is odd in [0,1] by Lemma 14.21 Since also sin 2 (^-i9 (s)) is even in 
[0, 1] by d4.4bl ). 04.341 * holds because its integrand is odd in [0, 1] . This concludes the proof 
of 04.28I *. Since the proof of 04.291 * is the same, we have shown that the second component 
of 04.261 * is zero. 

Therefore for every t e [0, to] there is a unique solution x\ (t) of 04.81 *. We now prove that 
x\(t) satisfies also 04.91 * for every t e [0,t ]- As before, this equation reduces to 04.13D . and 
again we will prove that both integrals in 04.131 * are zero for every t e [0, t Q ] . 

The proof for the first integral is the same as in the case t e [to, 1 — to] . Similarly, the 
velocity £(s, t) is zero if s ^ [L - 21, L] , so the second integral in 04.131 * reduces to an integral 
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on the inteval [L — 2l,L], After changing variables, we need to prove that 







L-21 



(K 6 (s,t)as,t),Jti(s,t))d S 



l 2 C T 

to 



cos(£0oOO) 
sin(^ (s)) 



tfo(<7)da 



(4.35) 



cos(£0 o (*)) 



sin(^ (s)) 



cos(£tf (<7)) 



da ) ds- 



l 2 C v 



cos(£tf (s)) 
sin(£tf ( S )) 



cos(£0 o (*)) 



da ) ds, 



which is satisfied if the following equalities hold true: 



cos 2 (^ (s)) 



sin(^-i?o(o-))da 



ds = 0, 



ds = 0, 



cos(£7?oOO)sin(£0 o OO) 



-l 
l 



S m(±Ma))M°)<la 



cos(^-^ (cr))d(T 



sin(^-^ (o-))dcr 



ds = 0, 
ds = 0, 



sin 2 (^ (s)) 



cos 2 (^ (s)) 



cos(^oW)i?o(<r)d<7 



cos(±Ma))M°)<lv 



cos{±ti (a)) da 



cos(^-i? (o-)) da 



ds = 0, 



ds = 0. 



This can be done as in the previous proofs, using the symmetry assumptions 04.4H together 
with Lemma l4T2l taking into account the intermediate identities 04.15D . d4. 17D . 04.22D . 04.28D . 

(Opt , and ma . 

Thus, equality 04.35D is proved, so that the solution ±i(i) obtained by solving 04. 8D also 
satisfies 04.9D in [0, t ] . 

The same proofs show that there exists a solution ±i(t) to 04.8D and 04.9& in [1 - t , 1] , 
therefore, we have obtained a unique solution ±i(t) to d4.8l > and 04. 9D in the whole interval 

[0,1]. 

4.1.3. Cyclic motion. We now estimate the overall translation x\(l) — xi(0) achieved by 
the maneuvers introduced in l4.1.1l and l4.1.2l We prove that it is strictly positive, provided 
the parameter I is suitably chosen. Then, in order to get and arbitrary prescribed transla- 
tion, we will repeat these maneuvers cyclically. 

First of all, we notice that the function xi(t) depends on the parameter I, related to the 
length of the bumps. We underline this dependence by using the notation x\ (t, I) . 

We begin with the estimate of the increment of xi(t, I) in the interval [to, I — h] , where 
the stationary regime takes place. The total viscous force Fl h (t) in the x± direction is (see 

(S3) 

(4.36) F° h (t) = -[ (Ks(s,t)Z(s,t),e 1 )ds, 

Jo 

and therefore, from 04.81 1 we obtain 



(4.37) 



xi (t,l) 



(K c (s,t)^(s,t), a) ds 
(K((s, i)ei,ei) ds 
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by taking into account d4.lt . 04. 2D . 04.7H . and 04.111 1 we get 

C T [ (l~cos§ Q (s))ds + (C u -C T ) I sin 2 ?9o(s)ds 
(4.38) ±i(t, I) = cl J —± 



la(l) 

sin 2 tf (s)ds 

for every t G [to, 1 - ^o] , where c is given in 04.5D and a(l) is a continuous positive function 
such that 

(4.39) lima(I) = a , 

i— >o 

for some a a e (0, +oo) . It follows that 

/•l-to 

(4.40) x 1 (l-t ,l)-x 1 (t ,l) = x 1 (t)dt = la(l). 

J to 

We now estimate the increment of xi(t, I) in the transient regime occurring in the inter- 
val [0, t ] . In this case, the total viscous force in 03.7H along x 1 is 

(4.41) *f(*)=- [ L (Ks(s,t)£(s,t),e 1 )ds, 

JL-21 

and thus 04.37D takes the form 



Xl(t,l) 



L 

(K i (s,t)i(s,t),e 1 )ds 

L-21 

L ' 

(K((s,t)ei,ei) ds 



recalling that £(s, t) = for s G [0, L - 21] . By JO ). (S3, <E7]>, and 04^261 ) we obtain 

z 2 c T f T (i) + a^(t) 



(4.42) Xl {tJ) 



tQ C T L + l(C v - C T ) I sin^ (s)ds 



where 



m - /; (/: ( ) m*)». ( %m )hnM^, 



Since F T and F„ are bounded, by 04.42H there exists a continuous positive function (3(1) and 
a constant /3 G (0, +oo) , with 

(4.43) Jim^(0=A), 
such that 

i±i(t,oi < l -m 

to 

for every te [0, to] . It follows that 

(4.44) |a;i(to,/)-a:i(0,0|</ 2 ^(0. 
Similarly, we can prove that 

(4.45) |a:i(l,0-a:i(l-to,0l < ^(0- 
By 04.4OH . 04.44D . and 04.45D we have 

(4.46) Axi(l) :=xi(l,i)-xi(0,0 G [ia(i) - 2J 3 0(J), la(i) + 2Z 2 /3(Z)]. 

By 04T391 . 04~43l >. and 04~46l ) there exists Z G (0, L/2) such that Axi(Z ) > 0. Moreover, 

(4.47) I \-> Ax-i(l) is continuous and lim Axi(l) = 0, 

by 0T381 . ( Q9l ). 04^21 . 0T43l >. and dOBl 

If the constant a > introduced after Lemma [472] satisfies a ^ Axi(7 ) , by 04.47D we can 
find a l a G (0, i ] such that Axi(Z a ) = a. Since we can always assume x±(0, 1) = 0, we obtain 
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that the equalities in ( 14.61 * are both satisfied for I = l a . This solves the controllability for 
the translation problem when a ^ Axi(Zq) • 

If a > Axi(l ), there exists an integer k ^ 1 such that a = kAxxQo) + a\ , where ^ 
ai < Axi(Iq). First, we construct the function £(s,i) corresponding to ^ on the rectangle 
[0, L] x [0, 1] and then we extend it by 1-periodicity in time to the rectangle [0, L] x [0, k] . The 
corresponding function \ will transfer the segment S in (s) = sei at time t = into the 
segment (kAxi(lo) + s)ei at time t = k. In order to transfer the latter segment into the 
segment £fi n (s) = (a + s)ei , in the time interval [k, k + 1] , we define £(s,t) on [0, L] x [fc, fc + 1] 
as in the previous step, with a replaced by ai . 

By the rate independence of the problems (see Remark 12.11 *. it is possible to rescale 
[0, k + 1] to [0, 1] . 

4.2. Rotation. In this subsection we describe how to rotate a straight swimmer about its 
center in the time interval [0,1]. This will be obtained in three steps. In the first one we de- 
form symmetrically the initial segment into the shape in Fig. [101 with two parallel straight 
terminal parts; by symmetry the deformation process will produce a rotation of an angle 
tpo (that we will not estimate) about the midpoint. In the second step we propagate bumps 
on the rectilinear parts as described below in order to achieve a rotation of a prescribed 
angle ip. In the third step, we straighten back the now rotated configuration in Fig.llOlinto 
a straight one by reverting the process in step one: this will produce a rotation of angle 
— tpo about the midpoint, so that at the end of the process the segment will be rotated by 
the angle ip . 




Figure 10. Rotation: steps 1 and 3 from the swimmer's reference system. 

Without loss of generality, in this section it is convenient to assume that the length of the 
swimmer is 2L and to parametrize all curves in the interval [— L, L] . We take Si n (s) = sei 
and Sfin(s) = se(if>), for s e \-L,L], where e((p) = (cos ip, sin ip) , ip being the angle of 
rotation. As before, the motion will be first described through a function £ £ X\ satisfying 
the two disks condition with radius p. Then, we will consider the function x(s,t) defined 
by 03.11 1. where x{t) and R(t) satisfy the equation of motion 03.8D . The initial and final 
conditions on x are 

(4.48) X {s,0) = se 1 , xO, 1) = se(tp), 

for all s £ [-L, L] . 

We also assume that £(s, 0) = sei and that £ satisfies 

(4.49) £(s,t) = -a~s,t) 

for every s e [-L, L] and for every t £ [0, 1] . It follows that 

(4.50) fat) = -i{-s,t), 

which implies that the force density f(s, t) in 02. ID is odd with respect to s, so that F(t) = 
$t L /(*. *) ds = for all t £ [0, 1] . 

The quantities introduced in 03.6D . 03. 7D . and 03.91 * are now defined by integration over 
the interval [— L, L] . The symmetry properties 04.49D and 04.501 1 imply also that the vector 
b(t) introduced in 03.6aD and the vector F sh (t) defined in 03.7D vanish. As a consequence, the 
vectors b(t) and v(t) introduced in 03.9H are zero and c(t) = l/c(t) . Therefore, the equations 
of motion 03. 8D read x(t) = and 

M sh (t) 

(4.51) 8(t) =u (t): = JE-m t 

c(t) 
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where 8(t) is the angle of the rotation R(t). Together with the initial conditions at time 
t = 0, this implies that x(t) = for every t G [0, 1] , 0(0) = 0, and 

X (s,t) = R(t)£(s,t). 

Therefore, the final condition in 04.48D is equivalent to 

0(1) = tp mod 2tt. 

The first step will take place in the time interval [0, 1/3] . The curve £(•, 1/3) is the one 
represented in Fig. [TO] The main feature of this curve, besides being odd, is that £(s, 1/3) , 
for s e [L/2,L], describes a segment on a straight line not passing through the origin. The 
angle ipo mentioned at the beginning of the sections is then defined by ip = 6(1/3) . 

The second step will take place in the time interval [1/3,2/3]. As in the case of pure 
translations, the overall rotation of angle tp will be achieved by iterating the cyclic motions 
described below. Since the problem is rate independent (see Remark I27IT ). each cycle will be 
parametrized in the time interval [0, 1] . 

During each cycle, we deform the rectilinear parts of the swimmer with the same bumps 
we used for the translation, see Fig. [9] They will be created in the time interval t <E [0, to] 
at the ends of the swimmer, will travel in the time interval t e [to, 1 — to] towards its center, 
and will be destroyed in the time interval t e [1 — to, 1] before entering the curvilinear part. 
As we proved in Section I4TT1 this motion will produce a net force F sh (t) directed along the 
rectilinear parts, see Fig. [TTJ 




P 

FIGURE 11. Rotation: step 2 (the right half of the swimmer is shown, s e [0, L\ ). 

As we obtained in 04.36D and in ( 14.411 1 the force F sh (t) generated during the translation 
of the bumps satisfies the estimates 

(4.52) F sh (t) = T^§-efa), for every t e [t , 1 - to], 

with a(l) -> a > as I -> 0, where e(rj) is the unit vector in the direction of the rectilinear 
parts, and 



f t0 \F sh (t)\dt+ f 1 \F sh (t)\dt^(3 l 2 , 

JO Jl-tn 



(4.53) 

-to 

for some f3 > 0. The corresponding torques defined in 03.7H can be computed as 

(4.54) M sh (t) = -2 [ L (J(Z(s,t)-p),K i (s,t)£( S ,t))d S + 2(jp,F sh (t))=2(jp,F sh (t)), 

JL/2 

where p := £(L/2, 1/3) is represented in Fig. [Til and the integral, which represents the 
torque with respect to the point p, vanishes, as explained in Section 1431 

To estimate 6 (t) in 04.511 1. we find bounds for c(t) from above and from below. Using 
fl3l3bl @3D, and (H~5Tfl we have 

c(t) = 2C T / (£'(s, t), J£(s, t)) 2 ds + 2C V \ {?(s, t), C(s, t)f ds. 

JL/2 JL/2 

By ( 13.101 *. 04. ID . and 04.3D . there exists constant < c min < c max such that 

(4.55) c min ^ c(t) ^ c max 
for every I small enough and for every tg [0,1]. 
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To estimate the angle 9(1) — 9(0) achieved in one cycle, we compute from 04.511 1 and 04.54D 

9(1) - 9(0) = 9(1) - 9(1 - t ) + 9(1 - t ) - 9(t ) + 9(t ) - 9(0) 

(4.56) f 1 M sh (t) , f 1 -* M sh (t) , f to M sh (t) , 

■ dt + / — T+ 1 dt + / A 1 dt. 



n-t c(t) A c(f) 7 c(t) 

First of all, we observe that, thanks to 04.531 1. 04.54D . and 04.55H . we can estimate 



(4.57) 



1 M sh (t) , f to M sh (t) , 
i-t c(t) Jo C W 



2 



To estimate the remaining integral in 04.561 1. it is not restrictive to assume that 

(Jp,e(rj))>0, 

with strict inequality since p and e(rj) are not collinear. From 04.521 1. 04.54D . and 04.551 1 we 
have 

(4.58) r° ^ dt = ^m^m r ta * dt > > . 

J t C (*) 1 - 2t O J t C (*) C max 

From 04361 1. 04371 ). and 04381 . we obtain that 9(1) - 9 (0) > , for Z small enough. 

Arguing as in the case of translations (see, Section [4. 1.31 1. we can prove that a rotation 
with prescribed angle <p can be obtained by iterating this procedure. 

At the end of step two the shape illustrated in Fig. [10] is rotated by an angle (p , there- 
fore, after scaling back step two on the time interval [1/3, 2/3] , we finally obtain 9(2/3) = 
9(1/3) + if = ifio + if. 

The third step will take place in the time interval [2/3,1]. We now define £(s,t) = 
£(s,l—t) for s e [—L,L] andfor t e [2/3,1]. Since this motion is the same as in step one with 
time reversed, the rotation angle 9(1) — 9(2/3) will be equal to — (po , hence 9(1) =<p. □ 

5. Optimal strokes 

In this section we prove Theorem 15.11 about the existence of an energetically optimal 
beating strategy. The result is achieved by proving that a minimum problem for the power 
expended 03.111 1 has a solution. 

Let us recall the definition of power expended: 

(5.1) V(x):= [ [ (-f(s,t),x(s,t)) dsdt= [ [ (K x ( S ,t)x(s,t),x(s,t)) dsdt. 

Jo Jo Jo Jo 

Up to a change of coordinates, it is possible to represent K x (s,t) in diagonal form, with 
entries C T and C v . Since < C T < C v , the matrix K x (s,t) is positive definite and its 
lower eigenvalue is C T . It follows that 

(5.2) V(x)>C T ( t \ X (s,t)\ 2 dsdt. 

Jo Jo 

For every p > let X{ be the set of all functions xe^i (see 02.12D ) such that for every 
t € [0, T] the curve x(-, t) satisfies the external disks condition with radius p (see Definition 

ED . 

Theorem 5.1. Let p > and let Xin, Xfm £ £T 2 (0,L;R 2 ), with \x' in (s)\ = \xL( s )\ = 1 f or 
every s e [0,L]. Assume that X'm and x&n satisfy the two disks condition with radius p (see 
Definition \3.4h Then the minimum problem 

(5.3) mm{V(x) : X £ X?, 023) holds, X (-,0) = Xin (-), x(;T) = Xfin(-)} 
has a solution. 

Proof. We first observe that the set of motions x ° n which we are mimimizing is nonempty 
by Theorem 14.11 Let us consider a minimizing sequence (xk)k- By 05. 2D there exists a 
constant M < +oo such that, 

(5.4) f [ \xk(s 7 1)\ 2 dsdt < M 

Jo Jo 
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for every k . Notice that the external disks condition with radius p gives the estimate 
(5.5) \xUs,t)\^l/p 

for every t G [0, T] and for a.e. s e [0, L] . 

We now show that Xfc(*>*) is bounded in L 2 (0,L) uniformly with respect to k and t. 
Since Xk(-, 0) = Xm(-)> f°r every s e [0, L] we have 

\ X k(s,t)\ 2 ^2|xi„( S )| 2 + 2T / |xfc(s,t)| 2 dt. 



From this inequality and from 05.4D we get 

Il 2 (0,L) * IIAin||£2( 0)i ) 



(5-6) sup ||xfc(-,*)ll^fo.w <2|| X in|li2 f0 , L) +2TM. 



te[o,T] 

By an elementary interpolation inequality we deduce from 05.5D and 05.6D that 

(5-7) sup \\Xk(;t)\\ H2{QX) ^C, 

te[o,T] 

for a suitable constant C < +oo independent of fc. 

By < I5.4I > and 05.61 1 the sequence (xk)k is bounded in i? 1 (0, T; L 2 (0, L)) . Therefore there 
exist a subsequence, not relabeled, and a function x '■ [0, £] x [0, T] — > M 2 such that 

(5.8) xefl" 1 (0 J T;Z a (0,L)), 

(5.9) Xk^X weakly in H\0 7 T; L 2 (0, L)). 

Since F^O, T; L 2 (0, L)) is continuously embedded into C°([0, T]; L 2 {0, L)) and for every t G 
[0,T] the function £ i-> £(t) is continuous from C°([0,T]; L 2 (0, L)) into L 2 (0,L), from dBT9l 
we deduce that 

Xfc(-,*) ^ *(-,*) weakly in i 2 (0,L) 
for every f G [0, T] . Then 05.7D gives that 

(5.10a) x(;t) 6 H 2 (0,L) for every i G [0,T], 

(5.10b) XfcM) xM) weakly in H 2 (0,L) for every t G [0,T], 

(5.10c) sup ||xO,i)l| ff 2 (OL) 

te[o,r] 

By (EB) and dSTTOcl) w e have x G X. Since the embedding of H 2 (0,L) into C 1 ([(),£]) 
is compact, from O5.10bl > we deduce that Xk(-,t) -> x(')*) strongly in C 1 ([0, L]) for every 
t 6 [0, T] . This allows us to pass to the limit in the equalities Xfc( s ' t) | = 1 , Xk(s, 0) = Xin(s) , 
Xk(s,T) = Xfin(s), and in the external disks condition with radius p. We conclude that 
X G X[ and that xO, 0) = Xm(s) , and x(s, T) = Xfln(s) . 

Let us verify that also the force and torque balance 02.3D passes to the limit. Equality 
02.3aH for \k reads 

(5.1D [ L K Xk (s,t)xk(s,t)ds = 0. 



Since Xfc('>*) converges to x'( - ^) strongly in C 1 ([0,i]) for every t G [0,T], by 02. 2D and 
05.1OcD we can apply the Dominated Convergence Theorem and we obtain 

(5.12) K Xk -> K x strongly in L 2 (0, T; L 2 (0, L)) . 
By 05.9H we have also 

(5.13) X k X weakly in L 2 (0, T; i 2 (0, i)) . 

By 05.11 * and by the Ioffe-Olech semicontinuity theorem (see, for instance, [6, Theorem 
2.3.1]) we have 

(5.14) V( X ) <liminfP(x fe ). 

k— »+oo 

Let now 95 G C c °(0, T) be a test function. By d5TTIt -( l5Tl3l we have 

<p(t)K Xk (s,t) X k(s,t)dsdt^ / / lf (t)K x (s,t)x(s,t)dsdt 



Jo Jo Jo 
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Since this equality holds for every ip e C°(0, T) , we conclude that 

/ K x (s,t) X (s,t)ds = 
Jo 

for almost every t e [0, T] . This proves 02.3aD . The proof of d2.3bP is analogous. Since (xk)k 
is a minimizing sequence, we deduce from 05.14D that \ is a minimizer of 05.3D . □ 
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